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We consider full counting statistics of spin in matrix product states. In particular, we study the 
approach to gaussian distribution for magnetization. We derive the asymptotic corrections to the 
central limit theorem for magnetization distribution for finite but large blocks in analogy to the 
Edgeworth series. We also show how central limit theorem like behavior is modified for certain 
states with topological characteristics such as the AKLT state. 



With the advantage of precision experiments in con- 
densed matter physics, it is now possible to probe the 
nature of correlated quantum systems to ever higher de- 
grees of detail and precision. This is especially apparent 
in cold atom systems, where in addition to refined cor- 
relation measurements, we also have fine control of the 
Hamiltonian parameters themselves. One of the most 
useful ways of studying detailed correlations in the in- 
termediate regime between the macroscopic, thermody- 
namic properties, and the microscopic, atom by atom 
level, is through the full counting statistics functions. 
These describe the full probability distribution of suit- 
able observables, such as the magnetization of a block of 
spins in a spin chain or the total excess charge flowing 
through a quantum point contact. 

The full counting statistics function (FCS) contains de- 
tailed information about the properties of the state. It 
has been a useful tool to analyze quantum states, from 
it's original appearance in quantum optics, in the theory 
of photon detectors [I] [5] in quantum optics to counting 
statistics of electrons in mesoscopic systems introduced 
by Levitov and Lesovik j3]. The full counting statistics 
has studied in numerous electronic systems theoretically 
[1HH] as well as in experiments [9J [TU] . The utility of full 
counting statistics for cold atoms was pointed out in [TT] . 
It has also been demonstrated that in certain cases count- 
ing statistics may be used to characterize block entangle- 
ment entropy in fermions states and spin states [T51U4) . 
Recently, the analyticity properties of the "bulk" compo- 
nent of the full counting statistics in classical Ising and 
quantum XY spin chains has been used as an alternative 
characterization of phases [15] . 

Here, we explore the quantum noise in an important 
class of Id states, the so called matrix product states 
(MPS) [TBJ. Such states appear as the ground states 
of frustration free Hamiltonians, and are especially use- 
ful as variational states to study other Id quantum sys- 
tems. MPS posses convenient properties that allow a 
thorough study of correlations and fluctuations in them, 
for example, an analogue of the Wick's theorem has been 
demonstrated for generic translationally invariant MPS 
in [TTJ. Finally, we remark that our results also hold for 
non quantum states as long as the probability distribu- 
tion of certain measurements may be described in terms 
analogous to MPS, a prominent example for which is the 



exact solution of the ID asymmetric exclusion process 
|18) . in a recent paper, the full counting statistics for the 
asymmetric exclusion model was considered in [19] . 

Here, we concentrate on the corrections to central limit 
theorem (CLT) like behavior of the full counting statis- 
tics. The central limit theorem is a description of the 
statistics of averages of independent random variables, 
stating that properly weighted average tend to a Gaus- 
sian distribution when the number of random variables is 
large. Since the correlation length in a MPS, say, is finite, 
one may expect gaussian like behavior for the magneti- 
zation of large blocks of spins. We find that the simple 
structure of the MPS allows us to not only do this but 
much more: we can controllably identify how the cen- 
tral limit of magnetization is reached, what are the main 
corrections (a consequence of entanglement in the sys- 
tem) and show how CLT may sometimes completely fail 
in cases of topological states. 

The distribution of magnetization approach to Gaus- 
sian behavior at large spin blocks is substantially more 
intricate for the MPS as opposed to independent random 
variables. To address this behavior we concentrate on 
deriving the asymptotic probability distribution and cor- 
rections to it. While in many cases, even when the distri- 
bution seems Gaussian in the infinite block size limit, the 
corrections due to finite block size are modified. Such cor- 
rections, are described, for independent, identically dis- 
tributed variables using various asymptotic series such 
as the Gram-Charlier A series and the Edgeworth se- 
ries pTJl HZT] . The Edgeworth series has been extensively 
studied in the mathematical literature, with the focus 
on ensuring it's applicability when dealing with random 
identically distributed variables, which may have diver- 
gent moments, see e. g. 

Here, we derive an expression for the asymptotic cor- 
rections to the central limit error function of MPS in 
analogy to the asymptotic Edgeworth series. We start 
by deriving formulas for the full counting statistics gen- 
crating function. For an extensive review of MPS their 
properties, and their relation to DMRG see [251 [SB]. MPS 
have been shown to be extremely well suitable for approx- 
imation of ground states of gapped Hamiltonian [27J HB] • 
Let us consider a MPS with periodic boundary conditions 
on N spins, defined as follows: 

i> = Z {a} Tr (A ai A a2 ..A aN ) \aia 2 -CT N ) (1) 



2 



where A a S {A%, ..As} , S is the spin index, and o\ e 
{1, ..S}. The matrices A are of size D x D, where D is 
called the bond dimension. 

To express the full counting statistics of the spin vari- 
able a, we define: 



E(X) = Y> a e lXa A a ® A 



(2) 



The full counting statistics generating function of the 
magnetization of a block of I sites is then given by 



x(A; I; N) = J2 n prob(total spin of block = n)e 



iXr 
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Here Si = Y,\ =1 <Ji where Oi is a spin operator at site i. 
When considering the thermodynamic limit, we add, as 
is usual, the demand that: 



^ a= \AaA^ — I 



(4) 



This ensures that the largest eigenvalue of E(0) is A = 1. 
In addition, when dealing with problems in the thermo- 
dynamic limit we assume this largest eigenvalue is non 
degenerate. In the thermodynamic limit, we define: 



X(A;Z)= lim x(A;Z;iV) 

N— too 



(5) 



The existence of the limit is assured by the conditions 
above. To compute x(A; I) we define -P(A) to be the ma- 
trix which brings E(X) to it's Jordan form, with Jordan 
blocks Jfc arranged such that the block with the largest 
eigenvalue is J\ . Note that the number of blocks as well 
as eigenvalues depend on A. We have: 

E{\) l =P{X){® k=l j{)p-\X) (6) 

We note that if there is no degeneracy, 

E(0) N ^ P(0)\1){1\P- 1 (0) (7) 

Therefore the full counting statistics function is given by: 

X (\;l) = (l\p- 1 (0)P(X)(® k=l J l k )p- 1 (X)P(0)\l) (8) 

Let a^(A) be the diagonal value of ■/&. We can compute 
explicitly the power of a Jordan block, obtaining after 
some algerba the formula: 

d k min(l,d k )—i 

X(A,0 = £E £ cl uQkM°k" ( 9 ) 

fc=l i=l u=0 

where C l v are the binomial coefficients, dk the dimension 
of Jordan block k and 



Qk,iAX) = {i + v + Yl* n =l d„|P- 1 (A)P(Q)|l) x 



(l\P- x {0)P{X)\i + v + Y!^ a d, 



(10) 



Note that in Q,ak as well as the limits in the sum 
depend on A implicitly. 

Since the largest eigenvalue 1 of E(0) is non- 
degeneracy, we have that Q\, 1,0(0) = 1 and Qk.i,v{fy = 
for all other value of k,i,v. It is also important note 
that since, generically, eigenvalues do not cross, we ex- 
pect that we may set dk = 1, v = in (J9J) for almost all 
values of A e [— tt,tt], unless some special symmetry or 
constraint is present. 

Let us now consider the limit of large block size 
As with any thermodynamic quantity, computed in 
a system with finite correlation length, we expect a 
gaussian distribution of observables according to the 
CLT. For a matrix product state, of course, the spins are 
not independent, and so the central limit distribution 
receives contributions from two types of corrections: due 
to correlations and due to finite size. Bellow we establish 
this behavior and derive the appropriate asymptotic 
description of the probability distribution for large but 
finite blocks. 

In the limit of large I, if the largest eigenvalue of the 
matrix E(X) is non degenerate, E(X) 1 is dominated by 
the largest eigenvalue ai(A) and we may write that: 

x(A;0~Xo(A)xi(A) i as I -> 00 (11) 



where: 



Xi(A) = ai(A) ; xo(A) = Qi,i,o(A) 



It is possible to take into account the corrections due to 
the smaller eigenvalues of E(X) as well, giving additional 
exponentially small corrections. 

We are now in position to describe the probability dis- 
tribution of block magnetization. Let us define: 



M, 



1 Si -Ifijl) 
yfl var(cr,l) 



(12) 



where fi(l) is the average magnetization per site, and 
var(a,l) is the variance per site. We note that since 
the spin variables on different sites are not independent, 
both (jl(1) and var(a, I) depend explicitly on the size of 
the block. Let 



Fi(M) = Prob(Mi < M) 



(13) 



be the probability distribution of measuring M{. To find 
F/(M), we now focus on the FCS for M;, defined as: 



Xm(A;0 



AXMi 



(14) 



Since we assume that the largest eigenvalue of E(X) is 
non degenerate at A — 0, this eigenvalue is analytic in a 
neighborhood of A = 0. Indeed both Xo(A), Xi(A) are an- 
alytic in the domain A < |A*| where A* is the smallest A 
(in the complex plain) for which the largest eigenvalue of 
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E(X) becomes degenerate (see, e.g. [H]). We can there- 
fore expand log(xi(X)), log(xo(X)) near A = 0. Noting 
that Xo(A) = Xl(A) = 1 we have the "cumulants" n r and 
£ r in: 



log(xo(X)) = £ 



r=l 

OO 



MXi(A)) = £ 



n! 

K r (iA) r 



(15) 



We see that for a block of I spins, fj,(l) — (a) — Ki + 
and var(a, I) = (a 2 ) — (a) 2 = ^/k 2 + £2//. We also 
recognize xo as boundary (or "Edge") term that char- 
acterize the effect of the rest of the chain on the chosen 
I spins, and X i as the bulk term that is not effected by 
other spins. 

In this paper xi pl a ys, formally, the role of the local 
independent random variable in the usual derivation of 
the central limit theorem. However, it is important to 
note that in a generic MPS, xi is n °t the full counting 
statistics of a valid probability distribution. Indeed, for 
that, the associated distribution, given by the Fourier 
transform of xi must be a positive real function. 

Let us briefly explore how close xi is to a valid proba- 
bility distribution. To do so we write the counting statis- 
tics of a block of size 1 as: 



X(A,1) =Xi(A)+X5(A), 
and define the pseudo-probabilities: 
1 



(16) 



Pn.i 



2tt 



dA Xl (A)e- 4An ; i = l,5. (17) 



Associated with the various eigenvalues of the 
Xi(A))X*(A). Thus p n< i is the effective probability 
distribution which would have generated the asymptotic 
behavior described by the central limit behavior. Wc 
can establish the following properties: 

1) From the definition of xi> we can immediately infer 
that the associated distribution is discrete. 

Indeed, observe that since E(X) is periodic, we can 
choose Xi(A) to be periodic: Xi(A) = Xi(A + 27r), which 
is associated with discrete, integer, spins. 

2) The distribution is real (however it is in general not 
necessarily positive). 

The second property is established by noting that: 



E(-X) = tE(X)t 
where r swaps Vi £3> V2 V2 ® v%. Therefore 



specE(-X) = specE(X) 



(18) 



(19) 



and, in particular Xi( — A) = Xi(A)j ensuring the Fourier 
transform is real. It is, however, in general not asso- 
ciated with a probability distribution, since the Fourier 
transform is, in general, not strictly positive. 



3 ) Enft.i = 1 and H n Pn,s = 0. To prove this note 
that xs(0) = since x(0, 1) = Xl(0) = 1. Now use that 

UnPn.S = X<5(0)- 

We now proceed to derive the Edgeworth series for our 



MPS. Using the definition ( 12 1 and cq. ( 11 1, we can find 
the cumulants for the distribution of Mi : 



logxM(X;l) 



(zA) 4 (Zk 4 + £ 4 



(iX) 2 (iA) 3 (fc 3 + &) 
2 6(Zk 2 +6) 3/2 

(iXf(lK 5 +£ 5 ) 



24(Zk 2 +£ 2 ) 2 120a« 2 +C 2 ) 5 / 2 "' 
We note that for the normal distribution, we have: 

(zA) 2 



(20) 



log(ct>(X)) = 



(21) 



Combine the two equations, and collect terms according 
to the power of I, we have 



log 



Xm(X) 1 (U) 3 k 3 , 1(zA) 4 k 4 



0(A) IV2 6k 3/2 

1 JiX) 5 K 5 (iX) 3 



p/2 



120k 



5/2 
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l 24k 2 
36 



2k 2 
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Exponentiate the above equation, we have 
Xm(A;Z) 



1 



\- gj(^A) , A 2 /2 



Z3/2 



(22) 



(23) 



where qj is a polynomial of degree 3j . 

Finally, to obtain FJ in (13), we do the inverse Fourier 



transformation to get the probability density, and inte- 
grate it over x to get the probability distribution. Defin- 
ing 



dq 

— ;=< 

2tt 



$(a;) = 

to be the error function. We obtain: 

Qj(-9 X ) 



3 = 1 



VI 2 



(24) 



(25) 



In general qj is a complicated polynomial, which can be 
compute to all orders. Here we write explicitly the first 
few terms: 



Q2 



6 k, 



3/2 



24k 2 



72k 3 , 



>4(d x ) 9 K 3 K 4 (d x ) 7 K 5 (d x f 



(26) 
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2k 2 ' 



4 



Comparing the above result with the usual Edgeworth 
series [20l |21] , we find the first two terms are the same as 
those appearing in the Edgeworth expansion for I inde- 
pendent measures with cumulants , the correction from 
the boundary term xo only effects the third and higher 
order terms. 

It is important to note that the parameters 
are, in principle, measurable. For example, K2,£2 can 
be obtained from the total noise in the measurement 
of I and I + 1 spins as k 2 = (A5 ; 2 +1 ) — (A6> ; 2 ) and 
6 = (I + l)(A5f) - l(ASf +} ), where ASj = 5, - (5/). 

Alternatively, by considering xoW as a differential op- 
erator acting on the Fourier transform of \\, and com- 
bining the Taylor series for xoW = Sfclo Ta^° with the 
Edgeworth series for x\i we may write explicitly 




Ft 



Prob( Sl 



< X 



Z^^Z^ K2 m/2 m | > 
L=Q m=0 z 



where Gk,m(%) is given by £7 ,m = (-d x ) m $(x): 



FIG. 1: Correction to the central limit distribution 
(i.e,F;(M) — <E>(M)). Solid line represent the exact result by 
calculating the probability distribution, dashed line shows the 
first order correction using Edgeworth series. 



Gk,m — 5] 

{pi, ••!>*} e z+ 

Sip; = fc; Sp; = j 



{-d x ) k+m+2j $(x) /K 3 \fl 



Pi!..p fe ! 



V3!/ "V( fc + 2 ) ! 



We can now compute FJ by using the expression (271, 
combined with: 



< x 



Prob [ Si=l£± < 



yf 



i+ 



i-2 



■Ft 



€2 V tK2 



To illustrate these ideas, let us consider the following 
spin 1 MPS, given by the properly normalized matrices: 



1/1 1 
3V0 



;A°= 



-1 1 
1 1 



1/0 
3V-1 -1 



Plugging these matrices in the definition ^ we find that: 



FIG. 2: First few Fourier components of xii showing small 
but finite negative pseudo-probabilities. 
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1 

-1 
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1 
1 



1 2e 



-iX 
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-iX 
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-iX 



V 



At A = we find that the largest eigenvalue is 1, and 
that it is separated by a gap from the next eigenvalue 
1/3. 

In Fig [T] we compare the probability distribution for 
magnetization computed numerically from the ground 
state wave function, with the probability distribution ob- 
tained from our Edgeworth series (25 1. Here, the exact 



probability distribution FJ (M) was found numerically by 
doing an inverse Fourier transformation of eq. ([3]). In 
the example depicted in Fig [l] for a block of 20 spins, 



it is evident that the Edgworth series works extremely 
well, capturing the essence of the correction already at 
first order. We also exhibit the pseudo-probabilities in 
Fig. [2j computed according to eq. (17 1. 

Next, we consider an example where the full counting 
statistics is not described by a Gaussian of finite width, 
although the system is gapped. In this example, the vari- 
ance per site actually vanishes as I /I. Consider the so 
called Affleck, Licb, Kennedy and Tasaki (AKLT) state. 
The AKLT state [30] is a unique ground state of the 
AKLT Hamilitonian. It has a "string", instead of local, 
order parameter, and also fractionalized edge excitations. 
See [3TJ 132] for a more recent study. The AKLT Hamil- 
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tonian is, 



^ = ^{S j S j+1 + ^(Sj + S j+1 ) 2 }. 



The ground state has an MPS form, according to 



A+ = 



2 /0 1 



3 00 



;A° = 



1 (-1 
3 I 1 




-1 



which gives: 



25(A) = 



1 



/ 1 



V2e 



2e iX \ 
-10 
0-10 
1 J 



Computing the eigenvalues of 25(A), we see that they 
are 1,— |,— g,— j, independent on A. In the limit of 
N,l oo, e s (A, Z, AT) LhpA, we find that the count- 
ing statistics does not depend on I. 

This result reflects the topological nature of the AKLT 
state, the total spin of the block depends only on the 
" edge modes" which are the only ones which are allowed 
to fluctuate. 

In general, the absence of scaling of fluctuations, ap- 
pears whenever we have a correspondence of the form 



25(A) = V(X)E(0)V- 1 (X) 



(28) 



For some matrix V. In such a case, Xi(A) = 1, and the 
entire contribution comes from the edge x(A; 1) = Xo(A). 
In this case logx(A) is clearly not extensive in the block 
size I. 

To summarize, Matrix product states supply a very 
natural class of probability distribution which are not 
IID, but are still quite amenable to treatment and com- 
putation. In this paper we have studied the full counting 
statistics of spin in matrix product states. We explored 
the finite size correction to the Gaussian counting statis- 
tics expected on large scales. We showed how an Edge- 
worth type series may be used describe these asymptotic 
corrections and checked it numerically on an explicit ex- 
ample. Finally, we showed that in special cases, such 
as the AKLT model, the fluctuations in the system do 
not scale linearly with system size, and the Edgeworth 
description is not valid whenever it is based on variance 
and mean which were measured on any finite size block. 
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